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ABSTRACT 

We present simulations of collapsing filaments studying the impact of turbulence and 
magnetic field morphologies on their evolution and star formation properties. We vary 
the mass per unit length of the filaments as well as the orientation of the magnetic 
field with respect to the major axis. We find that the filaments, which have no or a 
perpendicular magnetic field, typically reveal a smaller width than the universal width 
of 0.1 pc proposed by e.g. Arzoumanian et al. (2011). We show that this also holds in 
the presence of supersonic turbulence and that accretion driven turbulence is too weak 
to stabilize the filaments along their radial direction. On the other hand, we find that 
a magnetic field that is parallel to the major axis can stabilize the filament against 
radial collapse resulting in widths of 0.1 pc. Furthermore, depending on the filament 
mass and magnetic field configuration, gravitational collapse and fragmentation in 
filaments occurs either in an edge-on way, uniformly distributed across the entire 
length, or in a mixed way. In the presence of initially moderate density perturbations, 
a centralized collapse towards a common gravitational centre occurs. Our simulations 
can thus reproduce different modes of fragmentation observed recently in star forming 
filaments. Moreover, we find that turbulent motions influence the distance between 
individual fragments along the filament, which does not always match the results of a 
Jeans analysis. 
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1 INTRODUCTION 

The importance of interstellar filaments for star formation 
has been highlighted for the first time by Schneider & 
Elmegreen (1979). More recently, the advent of the Her- 
schel satellite has uncovered the filamentary structure in 
molecular clouds and infrared dark clouds in great detail 
(e.g. Andre et al. 2010; Konyves et al. 2010; Arzoumanian 
et al. 2011, 2013; Peretto et al. 2012; Schneider et al. 2012; 
Palmeirim et al. 2013, but see also the review of Andre et al. 
2014). 

The Herschel observations have revealed the interesting 
result that filaments apparently have a characteristic width 
of about 0.1 pc and a density profile, which can be repre¬ 
sented by a Plummer-like profile. Furthermore, line emission 
observations of e.g. C 18 0 or N 2 H + towards star forming fil¬ 
aments (Arzoumanian et al. 2013; Hacar et al. 2013; Furuya 
et al. 2014; Henshaw et al. 2014; Jimenez-Serra et al. 2014; 
Li et al. 2014) report non-thermal line broadening, which 
points to an additional trans- to supersonic turbulent ve¬ 
locity field. Temperatures in filaments are typically in the 
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range of ~ 10 to 15 K, which is particularly important for 
calculating the critical mass per unit length of a filament 
(Ostriker 1964). This is the mass per length below which 
the filament can be stabilized against radial (i.e. perpen¬ 
dicular to the major axis) gravitational collapse. Filaments 
with a mass per unit length above the critical value are grav¬ 
itationally unstable along the radial direction and will start 
to fragment along their major axis (e.g Myers 2009; Miet- 
tinen 2012; Kirk et al. 2013a,b; Zernickel et al. 2013). In¬ 
deed, observations have shown that filaments with observed 
star formation activity typically have a mass per unit length 
above the critical one, whereas filaments with a mass per 
unit length below this value are mainly devoid of any star 
formation (Andre et al. 2010). 

Polarisation measurements towards star forming fila¬ 
ments have also revealed the presence of large-scale and 
rather well-ordered magnetic fields. In many of the obser¬ 
vations the magnetic field appears to be roughly perpendic¬ 
ular to the filament (e.g. Chapman et al. 2011; Sugitani et al. 
2011; Planck Collaboration et al. 2014a). However, Li et al. 
(2013) report a bimodal distribution of magnetic field direc¬ 
tions (parallel and perpendicular orientation) with respect 
to the major axis of the filamentary structure (see also Pillai 
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et al. 2015; Planck Collaboration et al. 2014b). Moreover, by 
combining numerical simulations of magnetised, molecular 
clouds and synthetic polarisation maps, Soler et al. (2013) 
show that the relative orientation of the magnetic field also 
depends on the initial magnetisation of the cloud in which 
the filaments are forming. The knowledge about the rela¬ 
tive orientation would therefore allow to draw conclusions 
on the importance of magnetic fields during the formation 
of filaments: super-Alfvenic turbulence would cause strong 
compression, which results in magnetic fields parallel to the 
elongated structure (e.g. Padoan et al. 2001). In contrast, for 
a magnetic field guided gravitational contraction the field is 
preferentially perpendicular to the elongated structure (e.g. 
Nakamura & Li 2008). Furthermore, magnetic fields will also 
influence the subsequent evolution of filaments, a point we 
are particularly interested in in this work. 

The stability and fragmentation properties of filaments 
with respect to various aspects like e.g. turbulent motions, 
external pressure confinement, or accretion onto the filament 
have been studied analytically in a number of papers (e.g. 
Ostriker 1964; Inutsuka & Miyama 1992; Fischera & Mar¬ 
tin 2012; Pon et al. 2011, 2012; Toala et al. 2012; Heitsch 
2013). It has been found that magnetic fields generally have 
a positive effect on the stability of the filaments (Nagasawa 
1987; Fiege & Pudritz 2000; Heitsch 2013; Tomisaka 2014). 
However, these studies have to make several simplifications 
like e.g. considering an infinitely long filament or treating 
turbulence as an additional pressure term. 

There are also a number of numerical studies on the 
evolution and stability of filamentary structures. One of the 
first 2D numerical simulations were performed by Bastien 
(1983), which demonstrate some typical fragmentation prop¬ 
erties of filaments. Subsequently, more realistic 3D simula¬ 
tions showed the sensitivity of the fragmentation proper¬ 
ties of filaments to the initial conditions like the mass per 
unit length, or the presence of global rotation or magnetic 
fields (e.g. Bastien et al. 1991; Bonnell et al. 1992; Tomisaka 
1995). More recently, Burkert & Hartmann (2004) report 
the formation of prominent clumps forming at the edge of 
elongated, collapsing structures (see also Clarke & Whit¬ 
worth 2015). On molecular cloud scales, the formation of 
star forming filaments has been studied by means of tur¬ 
bulent box simulations and colliding flows (e.g. Mac Low & 
Klessen 2004; Hennebelle et al. 2008; Federrath et al. 2010b; 
Gomez & Vazquez-Semadeni 2014; Moeckel & Burkert 2014; 
Smith et al. 2014; Kirk et al. 2015). 

The present study systematically explores the impact 
of different magnetic field morphologies and turbulent mo¬ 
tions on the evolution of star forming filaments by means 
of 3D MHD simulations. We find that for supercritical fila¬ 
ments a magnetic field perpendicular to the major axis has 
almost no effect on its stability, which results in rather nar¬ 
row filaments. Only magnetic fields parallel to the major axis 
can stabilize the filament against radial collapse and reduce 
fragmentation. We will also show that turbulence strongly 
affects the fragmentation properties of filaments. 

The paper is structured as follows: in Section 2, we 
present an overview of the numerical methods and the ini¬ 
tial conditions used in the simulations. Next, we present the 
results focussing on the impact of the different initial con¬ 
ditions, before we investigate the fragmentation, the radial 
density profiles, and general star formation properties of the 


filaments. In Section 4 we discuss our findings and compare 
to recent observations and other numerical work, before we 
summarize in Section 5. 


2 NUMERICAL METHODS AND INITIAL 
CONDITIONS 

2.1 Numerical methods 

We present 3D, magnetohydrodynamical (MHD) simula¬ 
tions of magnetized filaments using the astrophysical code 
FLASH 4.2.2 (Fryxell et al. 2000; Dubey et al. 2008). We 
solve the equations of ideal MHD on an adaptive mesh 
including self-gravity. The Poisson equation for gravity is 
solved using a multipole method based on a Barnes-Hut 
tree 1 . The employed MHD-solver preserves positive states 
and is well suited for highly supersonic, astrophysical prob¬ 
lems (Bouchut et al. 2007; Waagan 2009; Bouchut et al. 
2010; Waagan et al. 2011). The maximum spatial resolu¬ 
tion, i.e. the size of the smallest grid cell is 40.3 AU. We 
apply a refinement criterion that guarantees that the Jeans 
length 



is resolved with at least 16 grid cells. Here c„ denotes the 
sound speed and G the gravitational constant. The gas is 
isothermal at a temperature of 15 K, which agrees well with 
temperatures observed in a number of filaments (e.g. Andre 
et al. 2010; Arzoumanian et al. 2011; Li et al. 2014). The 
mean molecular weight is set to /i lno i =2.3 typical for molec¬ 
ular gas. In order to guarantee that the Jeans length is re¬ 
solved everywhere, we create sink particles if the density 
exceeds a value of 

Psink = 10 _16 gcnU 3 (2) 

(for details on additional formation criteria see Federrath 
et al. 2010a). All gas in excess of that density, which is within 
a radius of 108 AU from the sink particle, is accreted. 


2.2 Initial conditions 


The properties of our filaments are motivated by number 
of recent observations (e.g. Andre et al. 2010; Arzouma¬ 
nian et al. 2011; Peretto et al. 2012; Busquet et al. 2013; 
Palmeirim et al. 2013; Schisano et al. 2014, but see also 
Andre et al. 2014 for a recent review). We set the initial 
filament width to 0.1 pc and apply a Plummer-like density 
profile along the radial direction: 


p(r) 


_Pc_ 

[1 + (J?/i?fl a t) 2 ] P/2 


( 3 ) 


Here R is the cylindrical radius, p c the central density, i?H a t 
the characteristic radius of the flat inner part of the density 
profile, and p the characteristic exponent of the profile. We 
note that the Plummer profile is applied for the entire radial 
range, in particular no cut-off density is applied. Arzouma¬ 
nian et al. (2011) find typical values for p between 1.3 and 
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2.4 with a mean of 1.6 and values for Ro a t between 0.01 pc 
and 0.08 pc with a mean of 0.03 pc, which gives a mean 
filament width of 3x ~ 0.1 pc. In this work we choose 
p — 2 and R^t = 0.033 pc. The length of the filaments is 
set to 1.6 pc. In order to avoid a pressure jump at the end 
of the filaments, the density decreases exponentially at each 
end of the filament. 

In Table 1 we list all simulations performed and their 
corresponding parameters. An important parameter for the 
evolution of a filament is its mass per unit length. As shown 
by Ostriker (1964), there is a critical mass per unit length 

(M/L )crit = ^ , (4) 

above which the filament starts to collapse along its radial 
direction and gets prone to gravitationally induced fragmen¬ 
tation ((A//L) C rit ~ 25 Mq/pc for gas with 15 K). In this 
work we investigate the evolution of filaments, which have 
( M/L)m = 1, 1.4, and 3 x (M/L) crit , i.e. ( M/L)m = 25, 35, 
and 75 Mq/pc (runs labelled with “FI”, “F2”, and “F3”, 
respectively, see Table 1). Given the density profile in equa¬ 
tion 3, this results in a central density of p c = 1 x 10 19 , 

1.4 x 1(T 19 , and 3 x 10 -19 g cm -3 , respectively. According 
to Arzoumanian et al. (2011), this can be converted into 
central column densities 


E c — A p p c .Rfl a t (5) 

with A p = ^ (1+ * 2 )p/2 ■ Assuming an inclination an¬ 

gle i = 0, this results in S c between 8.4 x 10 21 and 25.2 x 10 21 
end 2 , comparable to the authors findings. 

Observations show that filaments are most likely mag¬ 
netized with a direction of the magnetic field which is either 
perpendicular to the major axis 2 of the filament (e.g. Chap¬ 
man et al. 2011; Sugitani et al. 2011; Planck Collaboration 
et al. 2014a) or parallel to it (Li et al. 2013; Pillai et al. 2015; 
Planck Collaboration et al. 2014b). The preferred direction 
of the magnetic field in the densest parts, however, is still 
not entirely certain: Dust polarisation maps usually reveal 
a strong decrease of the polarisation towards the densest 
regions, which leads to some uncertainty about the field di¬ 
rection within the filament (e.g. Hildebrand et al. 1999; Hull 
et al. 2014). Furthermore, as shown by Soler et al. (2013), 
also the initial magnetisation of the surrounding molecular 
cloud can have a significant impact on the relative orien¬ 
tation of the magnetic field. Therefore, we will test both 
extreme cases here, i.e. perpendicular and parallel magnetic 
fields (run names contain “perp” or “para”; runs without a 
magnetic field are termed “NoMag”). We take the magnetic 
field strength to be 40 pG, which results in a magnetic pres¬ 
sure in the centre of the filament comparable to the thermal 
pressure, i.e. 


C s — VA , 

where va is the Alfven velocity 


B 



( 6 ) 

( 7 ) 


2 Throughout the paper we will use major axis for the long (sym¬ 
metry) axis of the filament and radial for the direction perpen¬ 
dicular to it. 


In the case of a perpendicular magnetic field this corre¬ 
sponds to a normalized mass-to-flux ratio between 1.6 and 
4.8 for the different simulations. Despite the fact that the 
choice of 40 pG is somewhat arbitrary, it compares reason¬ 
ably well with recent estimates of magnetic field strengths 
within filaments (e.g. Alves et al. 2008; Chapman et al. 2011; 
Sugitani et al. 2011). Moreover, as shown by Crutcher (2012, 
see their figures 6 and 7) at a volume particle density of 
n ~ 10 4 — 10 5 cm -3 or a column density around 10 22 cm~ 2 
as present in our filaments, the corresponding magnetic field 
strength in the galactic interstellar medium is of the order 
of a few 10 pG, which agrees well with our choice. The mag¬ 
netic field is initially constant throughout the simulation 
domain. Only for run Fl_perp_L2_Brad, the magnetic field 
decreases along the radial direction according to (Crutcher 
2012) 

B(R) = B 0 , (8) 

with the cylindrical radius R and Bo = 40 pG. 

Observations of star-forming filaments reveal non- 
thermal line widths, which indicate an additional turbulent 
velocity field with typical Mach numbers between 1 and 3 
(Arzoumanian et al. 2013; Hacar et al. 2013; Furuya et al. 
2014; Henshaw et al. 2014; Jimenez-Serra et al. 2014; Li 
et al. 2014). For this reason, we apply a transonic or super¬ 
sonic turbulent velocity field. We use a power-law spectrum 
with a slope of -11/3, corresponding to a Kolmogorov type 
spectrum in 3D. For each of the simulations, we choose a 
different random seed for the turbulence field. The typical 
length of the largest velocity fluctuations, i.e. the integral 
scale Amax of the turbulence field, is 0.1 pc, thus identical to 
the filament width. We change the turbulent Mach number 
and Amax, which is indicated by the run name (see Table 1), 
where “LI”, “L2”, and “L3” are used for A max = 0.0375,0.1, 
and 0.3 pc, respectively. In all simulations the smallest scale 
for velocity fluctuations is 80 AU, i.e. 2 grid cells, the Mach 
number (with respect to the sound speed) in all runs but 
one is 1. This is at the low side but still in agreement 
with the observational results mentioned above. We also per¬ 
formed an additional run with a Mach number of 2.5 (run 
F3_NoMag_L2_M2.5). We note that the Alfven Mach num¬ 
ber is comparable to the (sonic) Mach number (see Equa¬ 
tion 6). 

For all runs but one (run F3_NoMag_L2_condensed) the 
density is constant along the major axis of the filament. 
For run F3_NoMagJL2_condensed we enhance the density 
towards the centre of the filament by a factor of three: 

p(r,X) = PcX (l + (i? 2 /i? flat ) 2 ) + l + (R 2 + x*)/Rlj) ’ (9) 

where R is the cylindrical radius, x the distance along the 
major axis, p c = 3 x 1CP 19 g cm -3 , and the exponent p set 
to 2. In this run the magnetic field is zero, the turbulent 
Mach number is 1, and A max = 0.1 pc. 

In order to minimize boundary effects from the simula¬ 
tion domain, we place the filament in the centre of a domain 
which extents over 0.8 pc in the y- and 2 -direction and 2.4 
pc in the x-direction. Thus, the effective maximum resolu¬ 
tion is 4096 2 x 12288 cells. We use isolated boundaries for 
the gravitational potential and open (Dirichlet) boundary 
conditions for the gas. 
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Table 1. Performed simulations of filaments with a length of 1.6 pc and a width 0.1 pc. We list the run name, the central density p c , 
the mass per unit length (M/L )gp the magnetic field configuration (with respect to the major axis of the filament) and strength, the 
Mach number, and the turbulent integral scale A max - 


Run 

Pc [g cm 3 ] 

(M/L) m [M Q /pc] 

B-field 

B [pG] 

Mach number 

Vax [pc] 

Fl_NoMag_L2 

1 X 10“ 19 

25 

no field 

0 

i 

0.1 

Fl_para_L2 

1 x 10 -19 

25 

parallel 

40 

i 

0.1 

Fl_perp_L2 

1 x 10“ 19 

25 

perpendicular 

40 

i 

0.1 

Fl_perp_L2_Bweak 

1 x io- 19 

25 

perpendicular 

10 

i 

0.1 

F2_NoMag_L2 

1.4 x 10“ 19 

35 

no field 

0 

i 

0.1 

F2_perp_2 

1.4 x IO -19 

35 

perpendicular 

40 

i 

0.1 

F3_NoMag_L2 

3 x icr 19 

75 

no field 

0 

i 

0.1 

F3_para_L2 

3 x 10 -19 

75 

parallel 

40 

i 

0.1 

F3_perp_L2 

3 x icr 19 

75 

perpendicular 

40 

i 

0.1 

F3_NoMag_Ll 

3 x 10 -19 

75 

no field 

0 

i 

0.0375 

F3_NoMag_L3 

3 x 10“ 19 

75 

no field 

0 

i 

0.3 

F3_NoMag_L2_M2.5 

3 x 10 -19 

75 

no field 

0 

2.5 

0.1 

F3_NoMag_L2_condensed“ 

3 x 10“ 19 

75“ 

no field 

0 

1 

0.1 

F3_para_Ll 

3 x IO -19 

75 

parallel 

40 

1 

0.0375 

F3_para_L3 

3 x icr 19 

75 

parallel 

40 

1 

0.3 

F3_para_L2_Brad b 

3 x icr 19 

75 

parallel 

40 6 

1 

0.1 


“ central density enhanced by a factor of 3, see equation 9 
b radial dependence of the magnetic field, see equation 8 


3 RESULTS 

In this section we present the results of the simulations. First 
(Section 3.1), we focus on the time evolution of three fidu¬ 
cial runs Fl_para_L2, F3_para_L2, and F3_perp_L2. Next, 
we analyse the dependence of the fragmentation mode on 
global filament properties (Section 3.2) before we consider 
the evolution of the radial density profiles (Section 3.3). In 
Section 3.4, we investigate the distances between the frag¬ 
ments in the filaments, the global star formation rates are 
discussed in Section 3.5. 

We evolve each simulation for t ev0 i = 100 kyr - 250 
kyr after the first protostar has formed at ff orm . We follow 
each simulation over a long enough time to be able to draw 
reliable conclusions about properties like fragmentation and 
star formation rates. For this reason, t evo i is different for each 
simulation. Hence, also the final simulation times t en d = 
Iform + levoi change. For run Fl_perp_L2, there is no star 
formation at all within 2 Myr and we stop the run at this 
point. 

The formation time of the first protostar varies from 
95 kyr up to 1.4 Myr depending on the mass, the magnetic 
field configuration, and turbulence strength of the simulated 
filaments. Using the central gas density p c , the resulting free- 
fall times 

= \fSdc (10) 

are 210, 178, and 122 kyr for filaments with {M/L)d\ 
= 25, 35, and 75 Mg/pc, respectively (86 kyr for run 
F3_NoMag_L2_condensed). As can be seen from Table 2, 
where all time scales are listed, tf orm is larger than tfi by up 
to a factor of about 7 (run Fl_perp_L2J3weak). This is due 
to the fact that in particular for the filaments with (M/L)m 
= 25 Mg /pc the thermal pressure substantially counteracts 
gravity and free-fall is not a valid assumption any more. Also 


for the filaments with ( M/L)m = 75 Mg/pc, tf orm is larger 
than ffj by up to a factor of ~ 3, which we attribute to the 
combined effect of thermal and magnetic pressure. We note 
that using the mean density of the filaments within a cylin¬ 
drical radius of 0.2 pc, which is about a factor of 10 lower 
than p c , the corresponding ts is larger by a factor of ~ 3, 
which generally results in a better agreement between tform 
and tff. 

Due to the open boundaries, gas is allowed to flow into 
the simulation domain. The total mass inflow rates of the dif¬ 
ferent simulations range from 10 -6 to 10 -5 Mg yr -1 . Since 
these rates are about one order of magnitude smaller than 
the accretion rates onto the sink particles (see Section 3.5), 
we can neglect the impact of this additional mass increase 
(a few per cent) in the box. 


3.1 Time evolution 

First, we consider the time evolution of the runs Fl_para_L2, 
F3_para_L2, and F3_perp_L2, which demonstrate some basic 
differences that occur in the evolution of filaments with dif¬ 
ferent magnetic field morphologies and values of (M/L)hi. 
In Fig. 1 we show the column density of the three runs at 

t tform/2, tform; aild tend • 

For run Fl_para_L2 (left column of Fig. 1), the filament 
starts to contract along its major axis towards its geometri¬ 
cal centre whereas is remains stable along its radial direction 
(see also Section 3.3 for more details). At the outer edges 
of the filament overdense condensations form, which in turn 
collapse to form sink particles at tf orm = 1112 kyr. Over 
time, the sinks move towards the centre of the filament. Such 
a behaviour was previously found in numerical simulations 
by Burkert & Hartmann (2004) and was explained (semi-) 
analytically by a number of authors (Pon et al. 2011, 2012; 
Clarke & Whitworth 2015). Interestingly, ahead of the con- 
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Table 2. Free-fall time fgp calculated using the central density p c , formation time tf orm of the first protostar, simulation time t evo i elapsed 
after ff orm , total simulated time t em j, total time-averaged accretion rate, fragmentation mode, and number of fragments at f en( j. 


Run 

iff [kyr] 

£form [kyr] 

^-evoi [kyr] 

tend [kyr] 

Mstar [10 4 Mq yr 4 ] 

fragmentation mode 

# fragments 

Fl_NoMag_L2 

210 

1309 

250 

1559 

1.11 

edge-on 

3 

Fl_para_L2 

210 

1122 

250 

1373 

0.391 

edge-on 

2 

Fl_perp_L2 

210 

— 

— 

2000 

0 

— 

0 

Fl_perp_L2_Bweak 

210 

1398 

250 

1648 

0.602 

edge-on 

3 

F2_NoMag_L2 

178 

692 

200 

893 

1.70 

mixed 

18 

F2_perp_2 

178 

1239 

150 

1389 

0.43 

edge-on 

2 

F3_NoMag_L2 

122 

239 

150 

389 

5.01 

mixed 

24 

F3_para_L2 

122 

282 

150 

432 

1.30 

mixed 

7 

F3_perp_L2 

122 

264 

150 

414 

3.60 

uniform 

120 

F3_NoMag_Ll 

122 

240 

110 

350 

7.78 

mixed 

135 

F3_NoMag_L3 

122 

163 

200 

363 

2.78 

mixed 

17 

F3_NoMag_L2_M2.5 

122 

293 

100 

393 

6.13 

mixed 

28 

F3_NoMag_L2_condensed 

86 

96 

150 

246 

2.34 

centralised 

14 

F3_para_Ll 

122 

290 

150 

440 

1.12 

edge-on 

2 

F3_para_L3 

122 

179 

200 

379 

1.52 

mixed 

6 

F3_para_L2_Brad 

122 

284 

100 

384 

2.57 

mixed 

9 
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Figure 1. Gas column density and line-of-sight averaged velocity field (green arrows) at f = tf orm /2, tf orm , and t en< j (from top to bottom) 
for the runs Fl_para_L2, F3 para L2, and F3 perp 1/2 (from left to right). Black dots represent sink particles. The left and middle column 
demonstrates the impact of varying filament masses, whereas the middle and right column demonstrate the impact of different magnetic 
field morphologies. Note that the snapshots in each row are not taken at identical physical times but at comparable evolutionary stages 
(see Table 2 for the actual physical times). 


densation, i.e. at the side towards the centre, the gas is accel¬ 
erated outwards in the direction of the clump, in agreement 
with the findings of Clarke & Whitworth (2015). Consider¬ 
ing the inner parts of the filament, no further fragmentation 
occurs. Throughout the paper we will call this fragmentation 
behaviour, which is dominated by two fragments at each end 
of the filament, an edge-on collapse mode. 


Next, we consider the time evolution of run F3_para_L2, 
which has ( M/L)m = 75Mg/pc ~ 3(M/L) cr it (middle col¬ 
umn of Fig. 1). Also in this run, the column density profile 
does not appear to change significantly over time. However, 
due to the higher mass, the first sink particle forms signif¬ 
icantly earlier (ff orm = 282 kyr) than in run Fl_para_L2, 
although, also in this run, the two first particles form at 
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each end of the filament. However, more sink particles sub¬ 
sequently form along the major axis of the filament, which 
is why we call this collapse mode mixed. 

In the right column of Fig. 1 we show the time evolution 
of run F3_perp_L2, which has a magnetic field perpendicular 
to the major axis of the filament. Heavy fragmentation oc¬ 
curs along the filament with a significantly larger number of 
fragments compared to run F3_para_L2 (see Table 2). This 
is due to the lack of an additional magnetic pressure which 
stabilizes the filament against radial contraction in the first 
two simulations considered. Consequently, the filament in 
run F3_perpJL2 is much thinner and has a higher central 
column density. In run F3_perp_L2 the collapse does not 
proceed in an edge-on fashion. Fragmentation rather occurs 
randomly along the entire filament, with the first sink parti¬ 
cles forming closer to the centre. We call this fragmentation 
mode uniform. 

The different fragmentation modes also impact the dis¬ 
tribution of fragment masses. For run Fl_para_L2 the two 
sink particles have quite similar masses of 5.3 and 4.5 Mg. 
For run F3_para_L2 the two most massive particles are the 
outermost particles with masses of 7.4 and 5.9 Mg, whereas 
the remaining five particles have smaller masses between 0.9 
and 1.4 Mg. For run F3_perp_L2, the situation is quite dif¬ 
ferent. Here only one massive particle with 4.6 Mg exists. 
Furthermore, there are 10 more sinks with masses between 
1 and 3.3 Mg, 102 with masses between 0.1 and 1 Mg, 
and 7 light particles with masses below 0.1 Mg. We empha¬ 
size that in general runs with the same fragmentation mode 
(see Table 2) also reveal a similar distribution of protostellar 
masses. The actual distribution, however, might also depend 
on the spatial resolution of the simulation. In particular, we 
would expect a larger number of lighter particles with in¬ 
creasing resolution. However, we tentatively expect the gen¬ 
eral trend described above to hold even for a higher spatial 
resolution. 

Overall, we demonstrate the impact of different initial 
properties like the filament mass and magnetic field mor¬ 
phology on the fragmentation mode, the number of frag¬ 
ments, the filament width, and the distribution of fragment 
masses. We find that even small changes in the initial config¬ 
uration of the filament can cause characteristic differences. 
In the following sections we discuss these differences in more 
detail. 


3.2 Fragmentation mode 

In this section, we discuss the effect of the initial conditions 
on the fragmentation mode for all simulations with a fixed 
turbulent integral scale A ma x = 0.1 pc and a turbulent Mach 
number of 1 (see Table 1). Here, we only consider the final 
state (tend) since usually the time evolution of the individual 
runs is qualitatively similar to one of the three fiducial runs 
discussed in Section 3.1. We will show that filaments with 
(M/L)a i = (M/L ) CI it preferentially collapse in an edge-on 
fashion. For more massive filaments this is not the case, 
and only longitudinal magnetic fields help to stabilize the 
filament against heavy fragmentation. An overview showing 
all simulation discussed below at t = tend is given in Fig. 2. 


3.2.1 Filaments with ( M/L)si = (M/L ) CT it 

Similar to run Fl_para_L2, in simulation Fl_NoMag_L2 the 
formation of the sink particles initially occurs at the edges 
of the filament although a third sink particle forms about 
229 kyr after the first sink particle in the inner region, 
approximately 6100 AU offset from the centre. The over¬ 
all similar behaviour between the both runs concerning the 
fragmentation properties is not surprising since both have 
(M/L)fii = (M/L ) cr it and are thus marginally stable against 
fragmentation even without the presence of a magnetic field. 
Interestingly, the first sink particle in run Fl_para_L2 is 
formed somewhat earlier than in run Fl_NoMag_L2 ( about 
15% of fform, see Table 2), which could be due to the mag¬ 
netic field, which channels the gas along the longitudinal di¬ 
rection resulting in a somewhat faster accumulation of mass 
at the edges of the filament. 

In run Fl_perp_L2, where the magnetic field is oriented 
perpendicular to the filament, the formation of sink parti¬ 
cles is completely suppressed up to t en a = 2 Myr at which 
we stop the simulation. Here, the additional magnetic pres¬ 
sure stabilizes the marginally critical filament against grav¬ 
itational collapse, which is in good agreement with a mass- 
to-flux ratio of 1.6, i.e. a filament which is only marginally 
magnetically supercritical. However, in order to explore the 
stabilizing effect of the magnetic field in more detail, we re¬ 
peat this simulation with a four times weaker magnetic field 
of 10 /rG (run Fl_perp_L2_Bweak). As listed in Table 2, 
star formation now occurs at ff orm = 1398 kyr, i.e. approx¬ 
imately at the same time as in the runs Fl_NoMagJL2 and 
Fl_para_L2. Moreover, also for this run collapse occurs in 
an edge-on fashion with a third fragment forming 244 kyr 
after the first sink particle about 3200 AU away from its 
geometrical centre. 

3.2.2 Filaments with (M/L)m = 1.4 x (M/L) CI n 

Next, we analyse the fragmentation mode of the runs 
F2_NoMag_L2 and F2_perp_2 which have a mass per length 
of 35 Mg/pc ~ 1.4(M/L) crit . For run F2_NoMagJL2 the 
fragmentation mode is mixed where the fragments form first 
at the edges and subsequently along the entire filament. For 
run F2_perp_2, fragmentation is still dominated by an edge- 
on collapse and no further sink particles form along the fila¬ 
ment up to tend- It appears that even for marginally super¬ 
critical filaments ((M/L)m ^ 1.4 x (M/L) CI n) a magnetic 
field perpendicular to its major axis can stabilize it against 
heavy fragmentation or even completely suppress star for¬ 
mation as in run Fl_perp_L2. We note that in the outer re¬ 
gions of run Fl_perp_L2 the magnetic pressure has started 
to push material outwards causing the peculiar velocity field 
structure. 

3.2.3 Filaments with (M/L)a i = 3 x (M/L) CI i t 

In run F3_NoMag_L2, fragmentation occurs in a mixed 
mode, where the first sink particles form at either end of 
the filament. The number of fragments is significantly higher 
than in run F3_para_L2 (7 vs. 24 fragments) which we at¬ 
tribute to the stabilizing effect of the magnetic field in 
the latter case. Interestingly, compared to run F3_perp_L2, 
fragmentation seems to be be somewhat reduced in run 
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Figure 2. Gas column density and line-of-sight averaged velocity field (green arrows) at t = t en< j for the runs with A m ax = 0.1 pc and 
a turbulent Mach number of 1. Black dots represent sink particles. The figure demonstrates the impact of the mass per length of the 
filament and the orientation of the magnetic field on the fragmentation properties. 


F3_NoMag_L2. Comparing the kinetic energy content of 
F3_perp_L2 and F3_NoMagJL2 (not shown here), we find 
that in the former case turbulence decays somewhat more 
quickly in the beginning. Hence, we argue that the turbu¬ 
lent pressure counteracting gravity is smaller, which in turn 
results in a higher number of fragments in run F3_perp_L2. 

To summarize, increasing the mass per length of a fil¬ 
ament results in an enhanced fragmentation over the entire 
extent of the filament. A magnetic field parallel to the fil¬ 
ament helps to stabilize it against radial collapse and sub¬ 
sequent fragmentation. For perpendicular magnetic fields, 
only for marginally supercritical filaments (( M/L)m ^ 
1.4(M/L) cr i t ) the field can contribute to the stabilization 
of the filament. 


3.3 Density profiles 

Initially, the filaments have a radial density profile in agree¬ 
ment with recent observations (see Section 2). In the follow¬ 
ing we analyse how the radial profiles evolve and whether 
their functional form is retained. 

We display the radial density profiles for runs 
Fl_para_L2, F3_para_L2, and F3_perp_L2 in Fig. 3 at four 
different times (t = 0, tf orm /2, tf orm , and tend)- In order to 
reduce local fluctuations, we average the density along the 
major axis, excluding the edges, which have already frag¬ 
mented. For run Fl_para_L2 we use the inner 0.8 pc (i.e. 0.4 
pc in either direction from the centre) at t = tf orm /2, 0.4 pc 
at f = fform, and 0.2 pc at t = t en d; for the runs F3_para_L2 
and F3_perp_L2 we take the inner 1.2 pc, 1.0 pc, and 0.8 pc 
at t = tform/2, tform, and tend, respectively. 

For run Fl_para_L2 (left panel), the density profile re¬ 


mains similar to the initial profile. This is also valid for 
run F3_para_L2 (middle panel) despite a small increase in 
the central density. However, in run F3_perp_L2, the central 
density increases by about one order of magnitude. This 
indicates that here the thermal pressure of the gas is not 
sufficient to stabilize the filament against radial collapse, in 
agreement with Fig. 1. We note that the central density at 
tform (green line) is somewhat higher than at t en d since part 
of the gas is accreted onto the sink particles, whose mass 
we do not take into account when calculating the density 
profiles at later stages. 

We fit equation 3 to the density profiles shown in Fig. 3. 
The resulting values for P, and p c are listed in Ta¬ 

ble 3. In addition, we show the fit for t = fend (black dashed 
lines) in Fig. 3. For the runs Fl_para_L2 and F3_para_L2 
(left and middle panel), the profiles evolve only little and 
the fit parameters are comparable to the initial values and 
therefore agree reasonably well with those of Arzoumanian 
et al. (2011). For run F3_perp_L2 (right panel of Fig. 3), only 
the fit at t = tform/2 is still similar to the initial one, whereas 
we obtain a much more narrow profile at later times. 

We repeat this analysis for the final density profiles of 
several other runs (see Table 3) and plot p against Ra^t 
at tend in Fig. 4. Only filaments with an initially paral¬ 
lel magnetic field are compatible with observed values in¬ 
dicated by the grey shaded region. For the remaining runs 
the density profile becomes more narrow and/or declines 
more steeply. Also for more turbulent initial conditions (run 
F3JNoMag_L2_M2.5) a rather narrow filament develops, a 
problem we discuss in more detail in Section 4.1. Inter¬ 
estingly, for run Fl_perp_L2 the fitting values fall into the 
observationally determined range as well. However, unlike 
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Fl_para_L2 F3_para_L2 F3_perp_L2 



r / pc r / pc r / pc 


Figure 3. Averaged density profile along the radial direction of the filament for the runs Fl_para_L2, F3_para_L2, and F3_perp_L2 (from 
left to right). The times are identical to those in Fig. 1 (see also text). The dashed black lines represent the best fit (equation 3) at 
t = t en( j ( re <i lines). 


Table 3. Best fit parameter i?fl at , p, and p c from equation 3 for 
the radial density profiles of selected runs. 


Run 

Rfl at [PC] 

P 

p c [g cm 3 ] 

FIJSToMag_L2, t = t en( ± 

0.011 

10 

101 x 10” 19 

Fl_para_L2, t = tf orm /2 

0.036 

2.1 

1.02 x 10” 19 

Fl_para_L2, t = tf orm 

0.037 

2.0 

1.05 x 10” 19 

Fl_para_L2, t = t en ^ 

0.030 

1.7 

0.87 x 10” 19 

Fl_perp_L2, t = t end 

0.028 

2.1 

4.4 x 10” 19 

F2_NoMag_L2, t = t e n d 

0.0075 

5.6 

96 x 10” 19 

F2_perp_2, t = t end 

0.043 

3.4 

2.9 x 10” 19 

F3_NoMag_L2, t = t en( ^ 

0.015 

12 

45 x 10” 19 

F3_para_L2, t = tform/2 

0.021 

1.9 

6.0 x 10” 19 

F3_para_L2, t = tf orm 

0.016 

2.0 

11.2 x 10” 19 

F3_para_L2, t = t en( j 

0.017 

1.9 

8.2 x 10” 19 

F3_perp_L2, t = ff orm /2 

0.022 

2.0 

6.0 x 10” 19 

F3_perp_L2, t = tf orm 

0.0014 

1.8 

385 x 10” 19 

F3_perp_L2, t = t en d 

0.0041 

1.8 

33 x 10” 19 

F3JSToMag_L2_M2.5, t = t en d 

0.0090 

5.8 

86 x 10” 19 


the other filaments, here no stars are formed. We sup¬ 
pose that in this case the magnetic pressure supporting 
the filament against collapse in the longitudinal direction 
as well as the thermal support against radial contraction, 
i.e. (M/L)f{ i = (M/L) cr it, are just sufficient to maintain a 
filament width of about 0.1 pc. 

Overall, we find that whether or not the radial density 
profile of a filament fits the observed values partly depends 
on the considered time. At late times, all runs but those with 
a longitudinal magnetic field show a rather narrow filament 
width and a steep drop-off in disagreement with observa¬ 
tional results. 


3.4 Fragmentation spacing 

In the following, we analyse the mean spacing Ax between 
distinct regions of star formation, where distinct implies that 
we subsume close binary or multiple systems. We focus in 
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■ NoMag 
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- initial 
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0.01 
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Figure 4. Filament width (here represented with itfl a t) and ex¬ 
ponent p of the density profile at the end of the runs listed in 
Table 3. Red, green, and blue points correspond to filaments with 
an initial central density of 1 XlO” 19 , 1.4 XlO” 19 , and 3 XlO” 19 
g cm” 3 , respectively (labelled FI, F2, and F3 as in Table 1), 
squares, circles, and triangles to no magnetic field, a parallel, and 
a perpendicular magnetic field, respectively. The black asterisk 
represents the initial parameter setting. The grey-shaded rectan¬ 
gular gives the observational range found by Arzoumanian et al. 
(2011). Only the runs with a magnetic field parallel to the fila¬ 
ment axis are within that range. 


our analysis on the filaments with ( M/L)m = 3(M/L) C rit 
and compare A* with the Jeans length Aj (Equation 1), 
which for these runs is 0.94 pc. Here, we have assumed a gas 
temperature of 15 K and a gas density of 3 x 10 -19 g cm -3 . 
Furthermore, in order to take into account the stabilizing 
effect of a magnetic field, we can replace c 9 in equation 1 
by c 9 + v\, where va is the Alfven velocity (equation 7). 
Assuming that va — c s (see Section 2), we thus obtain a 
modified Jeans length, which is larger by a factor of V2. 

By means of a linear perturbation analysis, Nagasawa 
(1987) showed that the typical distance of fragments in 
a self-gravitating, isothermal cylinder with a longitudinal 
magnetic field is 

Acyl = 3.52 x (ttCs/Gpc) 172 . 
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Figure 5. Mean separation Ax between the sink particles for 
the runs F3_NoMag_Ll, F3_NoMag_L2, and F3_NoMag_L3 (left 
side), and F3_para_Ll, F3_para_L2, and F3_para_L3 (right side) 
after 350 kyr. The error bars show the minimum and maximum 
separation found in each run, the different colours and symbols 
distinguish the different turbulent integral scales A max . The two 
black dashed lines denote the Jeans length Aj excluding (lower 
line) and including (upper line) the magnetic pressure (see Sec¬ 
tion 3.4). The black line shows the spacing A cy i derived for an 
isothermal, self-gravitating cylinder. Since in run F3 para L! only 
two sink particles have formed, no minimum and maximum spac¬ 
ing is shown. 


Using T — 15 K and p c = 3 x 10 -19 g cm -3 , for these 
“sausage-type” instabilities we obtain a typical spacing of 
Acyi = 0.33 pc. 

We only consider the filaments with (M/L )fu = 
3(M/L) clit and determine the spacing 350 kyr after the start 
of the simulation. The results are shown in Figure 5 where, 
in addition to Ax, the error bar shows the maximum and 
minimum distance between the distinct sites of star forma¬ 
tion. For run F3_NoMag_L2, we identify 18 distinct locations 
of ongoing star formation. Given the fact that the distances 
shrink over time due to the global contraction of the fila¬ 
ment along its major axis, we normalize the spacing between 
the fragments to the initial length of the filament of 1.6 pc. 
Hence, for a given number of N star forming sites, we obtain 
the mean spacing by 


Ax 


1.6pc 
N - 1 ' 


( 12 ) 


For run F3_NoMag_L2 this results in a mean spacing of 
Ax ~ 0.094 pc. For run F3_perp_L2 (not shown in Fig¬ 
ure 5) we find about 15 different sites of star formation, 
which results in a mean spacing of about 0.11 pc, which is 
similar to that of run F3_NoMag_L2 despite the larger num¬ 
ber of fragments. In both runs Ax is in rough agreement 
with the initial Jeans length of the filament of Aj = 0.094 
pc. We point out, however, that this comparison is only ap¬ 
proximate since the density of the filament changes during 
its evolution. Nevertheless, this again demonstrates that a 
magnetic field perpendicular to the major axis of the fila¬ 
ment cannot contribute to the stabilization of a supercritical 
filament. 


For run F3_para_L2 the situation is somewhat different: 
Taking into account the overdense regions around x ~ —0.2 
pc and 0.05 pc (see middle panel in Fig. 1), which have not 
yet formed any sink particles, we obtain 8 regions of (future) 
star formation with a mean spacing Ax of about 0.23 pc. 
Considering only those regions which have actually formed 
sink particles, the mean distance Ax is about 0.32 pc. This 
is about 2-3 times larger than the distance of the frag¬ 
ments expected from a simple Jeans analysis, which gives 
Aj = 0.094 and 0.13 pc when excluding/including the mag¬ 
netic field pressure. The measured Ax, however, is in good 
agreement with the spacing of the “sausage-type” instability 
of A cy i = 0.33 pc (Nagasawa 1987). 


3-4-1 Impact of the integral scale of turbulence 

As shown before, the spacing of fragments can be explained 
by a Jeans analysis or the instabilities in an isothermal, self- 
gravitating cylinder (Nagasawa 1987). However, one could 
expect the turbulent motions within the filaments to also in¬ 
fluence the fragmentation properties, in particular since the 
turbulent integral scale A max = 0.1 pc is comparable to the 
mean distance between the fragments of 0.1 - 0.3 pc reported 
in Section 3.4. In order to test this, we have performed ad¬ 
ditional runs with turbulent integral scales A max of 0.0375 
pc and 0.3 pc both for no magnetic field and a longitudinal 
field (runs F3_NoMag_Ll, F3_NoMag_L3, F3_para_Ll, and 
F3_paraJL3) 

We plot the mean spacing Ax of these runs in Figure 5. 
For the unmagnetised runs there is a positive correlation 
between Ax and A max , which indicates a strong impact of 
the turbulent integral scale on Ax for unmagnetised, super¬ 
critical filaments 3 . Moreover, the poor agreement between 
Ax and the Jeans length for the runs F3_NoMag_Ll and 
F3_NoMag_L3 questions whether in the context of filament 
fragmentation the Jeans analysis is always a useful method. 

When excluding run F3_para_Ll, where only two sink 
particles have formed, for the magnetized runs there is no 
correlation between A max and Ax recognisable. The spacing 
is significantly larger than the Jeans length but in excellent 
agreement with A cy i (equation 11). Moreover, all values of 
Ax are larger than those in the unmagnetised runs, which 
again emphasizes the stabilizing effect of a longitudinal mag¬ 
netic field. 

A particular interesting point is the large difference in 
Ax for F3_NoMag_Ll and F3_para_Ll, which we attribute to 
the same physical origin: The smaller A ma x, the larger is the 
fraction of turbulent energy on smaller scales and the faster 
it gets dissipated (Stone et al. 1998). For run F3_NoMag_Ll 
this results in a quick decrease of the turbulent pressure 
counteracting gravity and in turn enhanced fragmentation. 
On the other hand, for run F3_para_Ll, where the mag¬ 
netic field already accounts for the stabilization of the fil¬ 
ament, turbulence appears to decay too quickly to be able 
to create local overdensities which could lead to additional 
fragmentation. Hence, the results of F3_NoMag_Ll and run 
F3_para_Ll represent the two modes of turbulence in star 


3 We note that due to the heavy fragmentation taking place in 
run F3_NoMag_Ll, it was hard to define a reliable minimum spac¬ 
ing (green arrow). 
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formation, i.e. globally counteracting gravity by means of an 
additional turbulent pressure and locally promoting star for¬ 
mation by creating local overdensities (Mac Low & Klessen 
2004). 

To summarize, our simulations show a high sensitiv¬ 
ity of the spacing between distinct sites of star formation on 
the turbulence field as well as the magnetic field orientation. 
The simple approach of a Jeans analysis does not match our 
findings for every case and could even by partly coinciden¬ 
tal. Rather, the filaments with a longitudinal magnetic field 
show a fragmentation behaviour in good agreement with the¬ 
oretical predictions for a self-gravitating, isothermal cylin¬ 
der. 


3.5 Effect on turbulence on star formation in 
filaments 

Next, we analyse the impact of the initial conditions on the 
global star formation properties. We show that our initial 
conditions have only a moderate impact on the global star 
formation rate although the previously discussed bimodal 
character of turbulence in star formation is recovered. 

In a first step, we consider the time averaged accre¬ 
tion rates M star of all sink particles and the formation times 
iform at which the first sink particles are formed (Table 2). 
Apart from run Fl_perp_L2 the values for M sta r and tf orm in 
each simulation-subset, i.e. runs with the same initial central 
density, do not differ too much from each other 4 . In fact, for 
each subset the total accretion rates lie within about half an 
order of magnitude and tform varies by less than a factor of 
two. 

The bimodal effect of turbulence on the star forma¬ 
tion process, i.e. locally promoting and globally hampering 
star formation (see e.g. Mac Low & Klessen 2004, for an 
overview, but see also Section 3.4.1), is also recognisable 
when considering the time when star formation sets in at 
tform: For the runs with (M/L )m = 3(M/L) cr i t , tform seems 
to decrease with increasing A max (see Table 2). We speculate 
that since for larger A max , turbulent velocity fluctuations 
occur also on larger scales, the turbulent energy does not 
dissipate as fast as in corresponding runs with smaller A max . 
The turbulent motions rather lead to strong local compres¬ 
sions of the gas which, in turn, results in a faster collapse 
and an earlier onset of star formation for larger A ma x- In 
contrast to that, when increasing the amount of initial tur¬ 
bulence from transonic to supersonic values (F3_NoMag_L2 
and F3_NoMag_L2_M2.5), it seems that the onset of star for¬ 
mation is somewhat delayed (tf orm = 239 kyr compared to 
tf orm = 293 kyr) due to a higher turbulent pressure. How¬ 
ever, the results also show that an enhanced amount of tur¬ 
bulence can by no means reduce the overall accretion rate 
M star once star formation has set in. Rather, there is a slight 
increase of about 20% in M s t ar with respect to M star of the 
fiducial run F3_NoMag_L2. However, since the difference is 
quite small, we cannot draw reliable conclusions. 


4 We exclude run F3_NoMag_L2_condensed from the considera¬ 
tion since here already the initial density profile differs. 


4 DISCUSSION 

4.1 A universal filament width? 

In the previous section we showed that for the runs without 
any magnetic field as well as a magnetic field perpendicular 
to the major axis of the filament, the fits of a Plummer-like 
profile at t = t en a fail to reproduce values comparable to 
observations (see Table 3 and Fig. 4), since the filaments 
become to too narrow. Only filaments with a longitudinal 
magnetic field reasonably fit the observations for all times. 
Furthermore, Fiege & Pudritz (2000) showed that also a 
toroidal magnetic field could stabilize the filament against 
radial collapse. However, observations show that the mag¬ 
netic field in filaments seems to be preferentially perpen¬ 
dicular to the major axis (Chapman et al. 2011; Sugitani 
et al. 2011; Planck Collaboration et al. 2014a) although also 
parallel configurations are observed (Li et al. 2013; Pillai 
et al. 2015; Planck Collaboration et al. 2014b). Under the 
assumption that there is indeed a universal width for star 
forming filaments of the order of 0.1 pc as proposed by Andre 
et al. (2014) (or even larger, see Schisano et al. 2014), our 
results raise the question what mechanism allows to retain 
this width. 

As suggested by Hacar et al. (2013), the apparently uni¬ 
versal filament width of about 0.1 pc could have its origin 
in the overlap of a number of parallel, rather narrow fila¬ 
ments, which are too close to each other to be resolved. By 
means of high-resolution, line emission observations the au¬ 
thors could show that filaments, which appear uniform in 
dust emission (e.g. Herschel observations by Arzoumanian 
et al. 2011), can split up into a bundle of narrow filaments 
with widths comparable to those found in our work. Such 
filamentary substructures were also seen in numerical simu¬ 
lations by Moeckel & Burkert (2014) and Smith et al. (2014). 
However, these authors follow the entire formation of the fil¬ 
aments in a turbulent environment, which probably sets the 
stage for an enhanced fragmentation process. In contrast, 
in our simulations we start from a pre-existing filament, 
which is most likely the reason why we do not find these 
sub-filaments. However, our filaments could represent one 
of these sub-filaments, indicating that such sub-structures 
would become rather narrow - as indeed observed. 


4-1.1 Accretion induced turbulence 

Arzoumanian et al. (2013) argue that during the evolution 
of the filament the turbulent support increases or is at least 
replenished due to accretion driven turbulence (Klessen & 
Hennebelle 2010; Goldbaum et al. 2011), which could sta¬ 
bilize the filament over time. This picture was confirmed 
by Heitsch (2013) by means of semi-analytical calculations. 
However, as demonstrated in Section 3.3, even for an ini¬ 
tially supersonic turbulence field (run F3_NoMag_L2_M2.5), 
the typical width of 0.1 pc cannot be retained for a super¬ 
critical filament, even though we allow mass to be accreted 
onto the filaments. 

In order to investigate this problem in more detail, we 
analyse the replenishment of turbulent energy in the case of 
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Figure 6. Mass infall rate (black line, left y-axis) and energy 
infall rate (red line, right y-axis) at the end of run F3 NoMag 1,2 
as a function of the cylindrical radius R. 


run F3_NoMag_L2 5 . In Fig. 6 we plot the mass infall rate 
Mfn(R) as a function of the cylindrical radius at the end 
of the simulation (black line). We obtain Mm{R) ^ 10~ 4 
Mg yW 1 . Furthermore, for the accretion driven turbulence 
model, the inflow rate of kinetic energy is a crucial quantity 
since it determines the amount of energy, which can - at 
least in parts - be converted into turbulent motions. The 
energy inflow rate is given by 


dE 

dt 


hi 


\Mm(R)v la . d (Rf , 


(13) 


where v ra d is the radial inflow velocity. We plot as a 

at | fil 

function of the radius R in Fig 6 (red line) obtaining values 
of the order of 10 30 — 10 32 erg s _1 . 

As shown by e.g. Mac Low et al. (1998) and Stone et al. 
(1998), the turbulent energy is roughly dissipated within one 
turbulent crossing time: 


dE_ 

dt 


diss 


^M(R)a 2 /tcioss = i M(R)a 2 /{R/a ), 


(14) 


where a is the mean turbulent velocity dispersion, R the 
radius which we set to 0.05 pc (half of the filament width 
of 0.1 pc), M(R) ~ 15.5 Mg the gas mass in the filament 
up to this radius, and r cross = R/a the turbulent crossing 
time. As stated in Klessen & Hennebelle (2010), only a small 
fraction e of the infalling energy can be converted into tur¬ 
bulent motions. Hence, in order to find out which amount 
of turbulence can be sustained by the observed infall rate, 
we have to solve the following equation for cr. 


dE 

dt 


1 , . <t 3 ! dE 

diss = 2 M(i?) It =e Ttt 


(15) 


Inserting the numeric values given above, we obtain a tur¬ 
bulent velocity dispersion of 


a ~ 0.5 km s 1 x e 1,/3 x 


dE I 
dt |fil 


10 31 ergs 



1/3 


(16) 


Assuming typical values for t between 10 3 and 10 1 and 
using the values for | from Fig. 6, we obtain values for a 


5 We note that the following estimate gives similar results for 
other runs as well. 


of at most 0.5 km s _1 , i.e. about 2.2 times the sound speed. 
Hence, with the observed infall rates, at most mildly su¬ 
personic turbulent motions can be sustained. This seems to 
agree with the results of Hacar et al. (2013) who show that 
(sub-)filaments in general appear to be “velocity-coherent” 
structures. However, with respect to a filament that is su¬ 
percritical, this turbulence level appears to be relative low 
and consequently in our simulations the filaments cannot be 
stabilized against radial collapse. 

A possibility to overcome this problem would be a more 
violent accretion flow, i.e. a larger Mm than the 10~ 5 to 10~ 4 
Mg yr“ 4 in the simulations presented here. This could be 
obtained by modelling a denser and gravitationally more un¬ 
stable molecular cloud environment as done by Smith et al. 
(2014) and Kirk et al. (2015). The authors find filament 
widths of the order of 0.1 pc and above, although Smith 
et al. (2014) report narrow substructures inside the filaments 
itself. 

To summarize, with the setup used here, supercritical 
filaments do not reveal a typical width of 0.1 pc at late 
stages unless a longitudinal magnetic field is present. An 
initial trans- to supersonic turbulent velocity field as well as 
the subsequent injection of turbulence by accretion (Klessen 
& Hennebelle 2010; Arzoumanian et al. 2013; Heitsch 2013) 
are not sufficient to stabilize a filament at a width of 0.1 
pc. To overcome this problem, turbulence would have to be 
replenished in a more efficient way, e.g. by higher inflow 
velocities than it is the case in our current simulations. Al¬ 
ternatively, filaments could be composed of a bundle of even 
more narrow, observationally unresolved (sub-)filaments or 
exhibit an additional longitudinal magnetic field component, 
which stabilizes the filament along its radial direction. 


4.2 Edge-on vs. centralised collapse 

Recently, Zernickel et al. (2013) observed a filament which 
exhibits two major, star forming clumps at either side. On 
the other hand, observations also show protostellar sources 
uniformly distributed over the entire filament (e.g. Arzou¬ 
manian et al. 2011; Miettinen 2012). In this work we showed 
that this could have its origin in different masses per unit 
length compared to the critical value (M/L ) CI it (equation 2). 
However, there are also a number of observations which 
report the collapse of filaments towards a central massive 
clump (Myers 2009; Schneider et al. 2012; Kirk et al. 2013a; 
Henshaw et al. 2014; Peretto et al. 2014), which we do not 
see in the simulations presented so far. A possible mechanism 
to obtain such a centralised collapse mode, is to impose dif¬ 
ferent initial conditions in our simulations. As demonstrated 
by Peretto et al. (2014) there seems to be an increase of the 
longitudinal, inwards directed velocity towards the common 
collapse centre. If there exists such an initially converging 
flow within the filament, this would explain both the velocity 
gradient as well as the mass accumulating at the convergence 
point of the flow. 

Another explanation would be that the collapse cen¬ 
tre, i.e. the massive clump in the centre of the filament 
was formed/present already during the formation process 
of the filament itself. Such a density enhancement could 
naturally lead to a collapse towards a common centre. 
In order to test this hypothesis, we consider simulation 
F3_NoMag_L2_condensed, where the density in the geomet- 
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Figure 7. Column density of the filament at the end of run 
F3_NoMag_L2_condensed. The increased central density results 
in a global collapse towards the centre. 

rical centre of the filament is enhanced by a factor of three 
(see equation 9). In Fig. 7 we show the column density of 
the filament at tend = 246 kyr. For this run, star formation 
already sets in at tf orm = 96 kyr, which is about three times 
earlier than in the corresponding run F3_NoMag_L2 (Ta¬ 
ble 2). Even for a moderate density enhancement of a factor 
of three, a centralised collapse occurs. Fragmentation in the 
outer parts occurs later, e.g. the sink particle at x ~ —0.4 
pc seen in Fig. 7 forms at t = 245 kyr. This agrees with 
the findings of Pon et al. (2011) (see their equation 24) who 
show that for a perturbation of size L\ and magnitude e 
in a cylinder of length 2 L, the collapse of the perturbation 
occurs faster than the collapse at the edge for 

L\ <C \/t x L . (17) 

Inserting Li = 0.1 pc, e = 2, and L — 0.8 pc as given in 
our simulation, the above condition is well met, in agreement 
with the observed centralized collapse mode. Since filaments 
naturally form in a turbulent environment, it is very likely 
that density fluctuations will occur along the major axis 
during their formation. Hence, a centralized collapse mode 
is easy to obtain for filaments which do not have a uniform 
density distribution along the major axis. 

4.3 Non-uniform magnetic fields 

So far the initial strength of the longitudinal magnetic field 
in our simulations was constant along the radial direction. 
In order to account for the observed scaling relation be¬ 
tween magnetic fields and the gas density (e.g. Crutcher 
2012) and to test the dependence of our results on the 
magnetic field configuration in more detail, we perform 
run F3_para_L2_Brad where the magnetic field strength de¬ 
creases along the radial direction (see equation 8). Compared 
to run F3_para_L2, the total accretion rate M s t ar on all sink 
particles is about a factor of 2 higher, which is caused by the 
reduced magnetic pressure. However, the creation time of 
the first sink particle as well as the number of sinks formed 
(7 in run F3_para_L2 and 8 in run F3_para_L2_Brad) are 
comparable, which leads to an overall similar fragmentation 
behaviour. 

Comparing the column density profile of 
F3_para_L2 J3rad to that of run F3_paraJL2 at t = 384 kyr 
shows that in run F3_para_L2_Brad the filament is more 
centrally condensed with a somewhat sharper outer edge 
(see Fig. 8). This is also confirmed by fitting the Plummer 
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Figure 8. Column density of the filaments of the runs F3_para_L2 
(top) and F3_para_L2_Brad (bottom) at t = 384 kyr. The initial 
radial decrease of the magnetic field in run F3_para_L2_Brad re¬ 
sults in a more condensed filament. 

profile (equation 3) to the radial density profile of the 
filament, which results in p c — 20 X 10 19 g cm 3 , i? flat 
= 0.014 pc, and p = 2.07. However, the values of Ra a t 
and p only differ by a few 10% from those of F3_para_L2 
and by a factor of 2 for p c . Since also the star formation 
properties agree reasonably well, the simplified assumption 
of a uniform longitudinal magnetic field seems sufficient. 

Although the initial conditions considered are highly 
idealised, the simulations can serve as a valuable guide for 
future and more sophisticated simulations including a better 
treatment of processes like chemistry and radiative transfer. 
We also point out that, due to the large number of simula¬ 
tions we had to perform for this study, an additional degree 
of freedom (some degree of disorder for example) would have 
not been feasible. Moreover, we emphasize that it is difficult 
to set up a more realistic filament, in particular when mag¬ 
netic fields are included, which requires the divergence-free 
condition to be fulfilled. In future work we will therefore in¬ 
vestigate the formation and evolution of filaments in molec¬ 
ular clouds in a self-similar fashion by using high resolution, 
zoom-in simulations of galactic discs (Walch et al. 2014). 
These simulations will also include a proper treatment of 
the thermodynamical properties of the gas, a fact we so far 
approximated by an isothermal equation of state. However, 
given the rather uniform temperatures found in interstellar 
filaments (e.g. Arzoumanian et al. 2011), we consider this as 
a reasonable approximation. 


5 CONCLUSIONS 

In this work we investigate the impact of turbulent motions 
and different magnetic field configurations on the evolution 
of interstellar filaments. We also vary the mass per unit 
length of the filaments ranging from one to three times the 
critical value. The simulations show a strong impact of tur¬ 
bulence and the magnetic field orientation on the evolution 
of interstellar filaments. They also provide a valuable guide 
for additional, more sophisticated simulations. 
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We consider the radial density profiles of the filaments 
and compare them with recent observations, which suggest 
a universal filament width of about 0.1 pc. In general we find 
that our simulated filaments have a much smaller width and 
a significantly steeper drop-off of the density at large radii 
than filaments observed in dust emission (see e.g. the review 
by Andre et al. 2014). In particular, magnetic fields perpen¬ 
dicular to the major axis cannot contribute to the stabiliza¬ 
tion of supercritical filaments ((M/L)m ^ 3(Af/L) cr it). We 
note that most of our filaments have rather moderate initial 
turbulence strengths (Mach numbers around 1), which are at 
the low side but still in agreement with the observational re¬ 
sults. However, even in the presence of supersonic turbulent 
motions, the filaments quickly become more narrow than ~ 
0.1 pc. We also show that accretion driven turbulence can 
only sustain at most mildly supersonic turbulent motions 
in our simulations. Hence, this mechanism is not sufficient 
to support the filaments along their radial direction and to 
retain a constant width unless a more efficient way is found 
to convert accretion energy into turbulent motions. In gen¬ 
eral, these findings fit into the idea that filaments could be 
composed of bundles of narrow (sub-)filaments (Hacar et al. 
2013). However, our simulations show that a typical width of 
0.1 pc can be retained if the magnetic field is oriented along 
the major axis of the filament which stabilizes it along the 
radial direction. 


Furthermore, we find that filaments with a mass per 
unit length equal to the critical one follow an edge-on col¬ 
lapse mode with star formation taking place at the outer 
edges of the filaments in agreement with numerical and 
(semi-)analytical results (Burkert & Hartmann 2004; Pon 
et al. 2011, 2012; Clarke & Whitworth 2015). No or only 
little fragmentation is found along the major axes of these 
filaments. Towards higher masses per unit length, more and 
more fragmentation is taking place along the entire filament 
(uniform collapse mode). When considering the evolution of 
a filament with an initial, moderate density enhancement in 
its centre (factor of three), we could show that it collapses 
towards this common gravitational centre ( centralized col¬ 
lapse mode). Our results are thus in good agreement with 
actual observations of star forming filaments, which discover 
all three collapse modes, edge-on, uniform, and centralized 
(e.g. Arzoumanian et al. 2011; Zernickel et al. 2013; Peretto 
et al. 2014). 


Finally, we find that the distance between fragments 
along the major axis can only partly be explained by a sim¬ 
ple Jeans stability analysis, even when the stabilizing effect 
of a longitudinal magnetic field is taken into account. For un¬ 
magnetised supercritical filaments we could show that tur¬ 
bulent motions have a strong impact on the fragmentation 
spacing. Hence, we tentatively suggest that the agreement 
between the fragmentation spacing and the Jeans length 
could be partly coincidental and that a simple Jeans analysis 
could be misleading. For filaments with a longitudinal mag¬ 
netic field, the fragmentation characteristics are reasonably 
well described by predictions for a self-gravitating cylinder 
(Nagasawa 1987). 
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